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1.  Introduction 

This  is  a  brief  report  on  the  work  done  by  the  IAU  working  group  ‘Relativity 
for  Celestial  Mechanics  and  Astrometry’  that  is  well  documented  on  the  Working 
Group  (WG)  web-site:  http://rcswww.urz.tu-dresden.de/~lohrmobs/iauwg.html. 
There  one  finds  a  list  of  members,  the  circulars  of  the  WG,  related  material  and 
references  to  relevant  publications  as  well  as  some  online  documents  containing 
important  formulas  and  explanations. 

The  first  central  task  of  the  WG  was  to  recommend  some  specific  form  of 
the  metric  tensor  gap  that  is  related  to  the  distance  ds  of  two  neighboring  points 
in  space-time  with  coordinates  xa  =  ( ct,xl )  and  xa  +  dxa  by 

ds2  =  gafj(t,  xl)  dxa  dx13 . 

The  metric  tensor  allows  one  to  derive  translational  and  rotational  equations  of 
motion  of  bodies,  to  describe  the  propagation  of  light,  the  rates  of  atomic  clocks 
and  to  model  the  processes  of  observation.  Meanwhile  it  is  widely  accepted  that 
in  order  to  describe  adequately  modern  astronomical  observations  one  has  to  use 
several  relativistic  reference  systems.  The  bary centric  celestial  reference  system 
(BCRS)  can  be  used  to  model  the  light  propagation  from  distant  celestial  objects 
as  well  as  the  motion  of  bodies  within  the  solar  system.  The  geocentric  celestial 
reference  system  (GCRS)  is  physically  adequate  to  describe  processes  occurring 
in  the  vicinity  of  the  Earth  (Earth’s  rotation,  motion  of  Earth’s  satellites). 

The  necessity  to  use  several  reference  systems  can  be  understood  from  the 
following.  If  we  were  to  characterize  terrestrial  observers  by  the  difference  be¬ 
tween  their  BCRS  coordinates  and  the  BCRS  coordinates  of  the  geocenter,  the 
positions  of  the  observers  relative  to  the  geocenter  would  change  with  time  also 
due  to  purely  relativistic  coordinate  effects  (such  as  Lorentz  contraction,  etc.) 
which  have  nothing  to  do  with  the  Earth’s  rotation  or  geophysical  factors  and 
vanish  if  one  employs  suitable  GCRS  coordinates  instead.  On  the  other  hand, 
the  coordinate  positions  derived  with  VLBI  observations  are  used  to  investi¬ 
gate  local  geophysical  processes  and  some  adequate  Geocentric  RS  allows  one 
to  simplify  their  description.  For  these  reasons  the  central  task  of  the  work¬ 
ing  group  is  to  specify  the  metric  tensors  both  in  the  BCRS  and  in  the  GCRS 
and  the  corresponding  space-time  coordinate  transformations  between  these  two 
systems. 

Two  advanced  relativistic  formalisms  have  been  elaborated  to  tackle  the 
problem  of  astronomical  reference  frames  in  the  first  post-Newtonian  approxi¬ 
mation  of  general  relativity.  One  formalism  is  due  to  Brumberg  and  Kopeikin 
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(Brumberg  and  Kopeikin,  1989;  Kopeikin,  1988;  Brumberg,  1991;  see  also  Klioner 
and  Voinov,  1993)  and  another  one  is  due  to  Damour,  Soffel  and  Xu  (Damour, 
Soffel,  Xu,  1991,  1992,  1993,  1994,  referred  to  as  DSX  I- IV).  Although  the  for¬ 
malisms  look  rather  different  at  first  glance,  it  is  mainly  the  concept  of  mass 
multipole  moments  (potential  coefficients)  at  the  first  post-Newtonian  level  that 
differs  in  the  two  formalisms.  The  new  recommendations  of  the  WG  improve 
and  extend  those  from  the  IAU  1991  framework  that  will  be  recalled  in  the  next 
section. 


2.  The  IAU  1991  framework 

The  IAU  resolution  A4(1991)  contains  nine  recommendations,  the  first  five  of 
which  are  directly  relevant  to  our  discussion. 

In  the  first  recommendation,  the  metric  tensor  in  space-time  coordinates 
(/,  x)  centered  at  the  barycenter  of  an  ensemble  of  masses  is  recommended  to  be 
written  in  the  form 


9oo  — 

-1  +  — 

got  = 

cs 

1 

CO 

9ij  = 

f  1  + 

+  0(c~% 


2U(t,  x) 


1  +  ^(c  4)  5 


(1) 


where  c  is  the  speed  of  light  in  vacuum  U  is  the  sum  of  the  gravitational  po¬ 
tentials  of  the  ensemble  of  masses,  and  of  a  tidal  potential  generated  by  bodies 
external  to  the  ensemble,  the  latter  vanishing  at  the  barycenter.  The  algebraic 
sign  of  U  is  taken  to  be  positive. 

This  recommendation  recognizes  that  space-time  cannot  be  described  by  a 
single  coordinate  system.  The  recommended  form  of  the  metric  tensor  can  be 
used  not  only  to  describe  the  barycentric  celestial  reference  system  of  the  whole 
solar  system  resulting  in  the  BCRS,  but  also  to  define  the  geocentric  celestial 
reference  system  (GCRS)  centered  at  the  center  of  mass  of  the  Earth  with  a 
suitable  function  U,  now  depending  upon  geocentric  coordinates.  In  analogy  to 
the  GCRS  a  corresponding  celestial  reference  system  can  be  constructed  for  any 
other  body  of  the  solar  system. 

In  the  second  recommendation,  the  origin  and  orientation  of  the  spatial 
coordinate  grids  for  the  solar  system  (BCRS)  and  for  the  Earth  (GCRS)  are 
defined.  Notably  it  is  specified  that  the  spatial  coordinates  of  these  systems 
should  show  no  global  rotation  with  respect  to  a  set  of  distant  extragalactic 
objects.  It  also  specifies  that  the  SI  (International  System  of  Units)  second  and 
the  SI  meter  should  be  the  physical  units  of  proper  time  and  proper  length  in 
all  coordinate  systems.  It  states  in  addition  that  the  time  coordinates  should  be 
derived  from  an  Earth  atomic  time  scale. 

The  third  recommendation  defines  TCB  (Barycentric  Coordinate  Time) 
and  TCG  (Geocentric  Coordinate  Time)  —  the  time  coordinates  of  the  BCRS 
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and  OCRS,  respectively.  Here  we  write  (t  =  TCB,xv)  and  (T  =  TCG,Xr) 
for  the  respective  coordinates.  The  recommendation  also  defines  the  origin 
of  the  time  scales  (their  reading  on  1977  January  1,  0h0m0s  TAI  ( .JD  = 
2443144.5  TAI)  must  be  1977  January  1,  0h0m  32.184s)  and  declares  that  the 
units  of  measurements  of  the  coordinate  times  of  all  reference  systems  must  co¬ 
incide  with  the  SI  second  and  SI  meter.  The  relationship  between  TCB  and 
TCG  is  given  by  a  full  4-dimensional  transformation 


TCB  -  TCG  =  c~2 


t  I  2  +  )  dt 


vErE 


0(c 


— 4> 


(2) 


where  xlE  and  vlE  are  the  barycentric  coordinate  position  and  velocity  of  the 
geocenter,  rE  =  xt  —  xE  with  xl  the  barycentric  position  of  the  observer,  and 
U(t,XE(t))  is  the  Newtonian  potential  of  all  solar  system  bodies  apart  from  the 
Earth  evaluated  at  the  geocenter. 

In  the  fourth  recommendation  another  time  coordinate,  Terrestrial  Time 
(TT),  is  defined  for  the  GRS.  It  differs  from  TCG  by  a  constant  rate  only 


TCG  -TT  =  Lgx  (JD  -  2443144.5)  X  86400,  LG  «  6.969291  X  10“10,  (3) 

so  that  the  unit  of  measurement  of  TT  agrees  with  the  SI  second  on  the  geoid. 
TT  represents  an  ideal  form  of  TAI ,  the  divergence  between  them  being  a  con¬ 
sequence  of  the  physical  defects  of  atomic  clocks.  It  is  also  recognized  that  the 
TT  is  nothing  else  than  a  rescaling  of  the  GRS  coordinate  time  TCG. 

The  fifth  recommendation  states  that  the  old  barycentric  time  TDB  may 
still  be  used  where  discontinuity  with  previous  work  is  deemed  to  be  undesirable. 
In  the  notes  to  the  third  recommendation  the  relation  of  the  TCB  with  TDB 
given  as 


TCB-TDB  =  Lbx(JD- 2443144.5)  X  86400,  Lb  ~  1.550505  X  10-8.  (4) 

3.  General  framework  for  new  conventions 

3.1.  New  conventions  for  the  Barycentric  Celestial  Reference  System 

The  metric  tensor  of  the  Barycentric  Celestial  Reference  System  (BCRS)  in  the 
first  post-Newtonian  approximation  should  be  written  in  the  form 


ffoo 

9ot 

9ij 


1 


(5) 
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Here,  the  post-Newtonian  gravitational  potential  w  generalizes  the  usual  New¬ 
tonian  potential  U  and  t»'  is  the  vector  potential  related  with  gravito-magnetic 
type  effects.  The  best  reason  for  writing  the  metric  in  that  form  is  a  simplic¬ 
ity  argument.  This  is  the  most  compact  form  to  write  the  metric  tensor  to 
first  post-Newtonian  order.  Let  us  also  note  here  that  the  potential  w  contains 
also  explicit  post-Newtonian  terms  and,  therefore,  does  not  coincide  with  the 
Newtonian  potential.  As  is  well  known  it  makes  no  sense  to  distinguish  be¬ 
tween  the  Newtonian  potential  and  explicit  post-Newtonian  terms.  It  is  the 
post-Newtonian  potential  w  as  a  whole  that  plays  a  role  in  observations. 

Note  that  this  form  (5)  of  the  barycentric  metric  tensor  implies  the  barycen- 
tric  spatial  coordinates  ad  to  satisfy  the  harmonic  gauge  condition.  We  rec¬ 
ommend  also  to  use  the  harmonic  gauge  for  the  barycentric  coordinate  time 
t  =  TCB.  The  main  arguments  in  favour  of  the  harmonic  gauge  are: 

•  tremendous  work  on  General  Relativity  has  been  done  with  the  harmonic 
gauge  that  was  found  to  be  a  useful  and  simplifying  gauge  for  all  kinds  of 
applications; 

•  in  contrast  to  the  standard  PN-gauge  the  harmonic  gauge  can  be  defined 
to  higher  PN-orders,  and  in  fact  for  the  exact  Einstein  theory  of  gravity. 

Assuming  space-time  to  be  asymptotically  flat  (no  gravitational  fields  far 
from  the  system),  i.e., 


lim  gtlv  =  diag(-l,+l,+l,+l) 

t  =  const 

in  the  standard  harmonic  gauge,  the  post-Newtonian  held  equations  of  General 
Relativity  are  solved  by 


w(t,x)  =  G  J  ffV  |^’^|  +  J  d3x'  a(t,x')\x  -  x'|  , 

w\t,  x)  =  G  j  d3x'^  X/| ,  (6) 

where 

rpOO  _i_  rpss  rpOi 

<j(t,x)=  - - - ,  ff'(fx)=  -  . 

cz  c 

T^v  =  Tlll/(t,xl)  are  the  components  of  the  energy-momentum  tensor  in  the 
barycentric  coordinate  system  and  Tss  =  T11  -\~T22  +T33.  For  many  applications 
explicit  expressions  for  the  gravitational  mass  and  mass-current  density,  a  and 
<7®,  will  not  be  needed. 

Since  the  integrations  in  eqs.  (6)  have  to  be  taken  over  all  the  massive  solar 
system  bodies  A,  the  metric  potentials,  w  and  w\  can  be  written  as  a  sum  of 
the  form 


w  =  J2wa'^ 

A 


w%  =  WA  ■ 
A 
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Here  the  potentials  with  index  A  are  obtained  from  relations  (6)  with  integrals 
taken  over  the  support  of  body  A  only.  In  the  following  we  will  use 

Wext  =  W Wext  =  WA  > 

AgE  AgE 

where  E  stands  for  the  Earth. 

3.2.  New  conventions  for  the  Geocentric  Celestial  Reference  System 

The  metric  tensor  of  the  Geocentric  Celestial  Reference  System  (GCRS)  (T,  X1) 
should  be  written  in  the  same  form  as  that  of  the  BCRS: 


G  oo 
Goa 

Gab 


2W  2 IT2  _5, 
-!  +  -2 - ZT  +  °(c 


4 

~~Wa, 

Sab  (l  +  +  0(C~4)  . 


(7) 


Here  IT  =  IT(T,  X)  is  the  post-Newtonian  gravitational  potential  in  the  geo¬ 
centric  system  and  IT“(T,  X)  is  the  corresponding  vector  potential.  These  geo¬ 
centric  potentials  should  be  split  into  two  parts:  potentials  We  and  WE  arising 
from  the  gravitational  action  of  the  Earth  and  external  parts  lText  and  IText  due 
to  tidal  and  inertial  effects.  The  external  parts  are  assumed  to  vanish  at  the 
geocenter  and  admit  an  expansion  into  positive  powers  of  X.  Explicitly, 


W (T,  X)  =  We(T,  X)  +  Wext(T,  X) , 
1T“(T,  X)  =  WE(T,  X)  +  lT“xt(T,  X) . 


(8) 


The  Earth’s  potentials  We  and  WE  are  defined  in  the  same  way  as  we  and  wlE 
but  with  quantities  calculated  in  the  GCRS.  We  may  write 


khext  —  Winer  +  ITtidal  7 


IT 


ext 


=  !T“  +  IT 

*  y  mer  '  r  y  t 


tidal  ' 


(9) 


Here,  lTiner  and  Wj“er  are  inertial  contributions  that  are  linear  in  Xa 
fCiner  =  Qa*\ 

=  lc2Sabc(nb-nler)XC+0(c-2). 


(10) 


Qa  characterizes  the  deviation  of  the  actual  worldline  of  the  origin  of  the  GCRS 
from  geodesic  motion  in  the  external  gravitational  held  that  is  determined  mainly 
by  the  coupling  of  the  Earth’s  nonsphericity  to  the  external  potential.  To  New¬ 
tonian  order  Qa  is  given  by 


Qa  =  Ra 


_d_ 

dxl 


wext  -atE)+  0(c-2) 


(11) 
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Here,  xE(t),vE(t)  =  dxE/dt  and  aE  =  dvE/dt  are  the  barycentric  coordinate 
position,  velocity  and  acceleration  of  the  origin  of  the  GCRS  (geocenter).  f?“ 
is  a  rotation  matrix  that  determines  the  relative  orientation  of  barycentric  and 
geocentric  spatial  coordinate  lines.  0“  is  a  slowly  varying  function  of  time  related 
to  f?“  by 


eabcttb  =  RakR%  ,  (12) 

and  fi“ner  mainly  describes  the  geodetic  precession  of  inertial  axes  with  respect 
to  remote  objects.  One  sees  that  for  0“  =  fi“ner  the  vector  potential  Wj“er 
vanishes.  This  implies  that  dynamical  equations  of  motion,  e.g.,  for  a  satellite 
around  the  Earth,  do  not  contain  Coriolis  and  centrifugal  terms,  he.,  the  local 
geocentric  spatial  coordinates  Xa  are  dynamically  non-rotating.  Recommended, 
however,  is  the  use  of  kinematically  non-rotating  geocentric  coordinates  defined 
by  tta  =  0,  be.,  by  Rat  =  6ai. 

pptidal  -g  a  generalization  of  the  Newtonian  tidal  potential 

X)  =  Wext(*E  +  X)  -  Wext(*E)  ~  Xa^LWext(*E)  •  (13) 

Full  post-Newtonian  expressions  for  ITtidal  and  J^tidal  can  found  in  DSX  II, 
IV. 

Finally,  the  local  gravitational  potentials  WE  and  WE  of  the  Earth  are 
related  to  the  barycentric  gravitational  potentials  we  and  wE  by 

We(T,  X)  =  wE(t,  x)  +  ^  ^  vlE  wlE(t,  x)  +  0(c~4), 

We(T,  X)  =  R?  (wUt,x)-vhwE(t,x))  +  0(c-2).  (14) 

3.3.  Transformations  between  the  reference  systems 

The  coordinate  transformations  between  the  BCRS  and  GCRS  can  be  written 
as 


T 


Xa 


t 


'  ~2  +  VErE ) 

+  1  (B(t)  +  B\t)  rE  +  Bij  (t)  rE  r3E  +  C(t,  x) 
+0(c~ 5), 

rE  +  ^2  (\vEvErE  +  WeA*E)rlE  +  r'Ea3Er3E  - 


+  0(c~% 
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:aErE 


(15) 


(16) 


2VE  +  Wext  (X£j), 


where 
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-\VE~\VE  Wext  (xB)  +  4  VlE  Wlext 


B\t)  =  -  i  v%  vlE  +  4  <xt(x£;)  -  3  vlE  wex t(xE), 

Bl3{t )  =  -v'E  Raj  Qa  +  2 ^J<xt(xE)  -  v'e  ^jWext(xs)  +  ^  <W  Wext(xS), 

<^(^x)  =  --ffirEiahrh)- 

Let  us  remark  that  the  harmonic  gauge  condition  does  not  fix  the  function  C 
uniquely.  Here  we  have  indicated  the  simplest  solution.  Though  the  theoretical 
domain  of  validity  of  the  GCRS  coordinates  is  quite  large  (larger  than  the  size 
of  the  solar  system)  the  spatial  region  where  the  GCRS  coordinates  are  to  be 
used  in  practice,  should  be  restricted  to  the  immediate  vicinity  of  the  Earth.  We 
suggest  here  that  the  GCRS  coordinates  be  used  only  to  about  the  geostationary 
orbit,  that  is,  for  |X|  <  50000  km.  The  size  of  this  region  allows  one  to  estimate 
the  terms  in  the  transformations  and  to  neglect  those  terms  which  are  smaller 
than  the  targeted  accuracy.  Thus,  it  is  easy  to  see  that  in  the  transformations 
between  coordinate  times  c~4  C  <  0.1  ps  for  |X|  <  0.1  AU  and  can  be  neglected 
for  most  applications.  In  the  same  way,  one  can  show  that  c~4BlJ  rErE  <  0.1 
ps  for  | X |  <  0.01  AU  and  also  can  be  neglected  for  most  purposes.  However,  to 
avoid  ambiguities  one  should  remember  that  the  terms  0(rE)  are  fixed  in  the 
time  transfromations  and  neglected  only  because  of  their  small  numerical  values 
in  the  considered  region  of  space. 

3.4.  Multipole  expansions  of  the  local  gravitational  potentials 

For  many  problems  it  is  advantageous  to  present  the  local  gravitational  poten¬ 
tials  of  the  Earth  as  multipole  series  that  usually  converge  everywhere  outside 
the  Earth.  The  definition  of  corresponding  post-Newtonian  multipole  moments 
or  potential  coefficients  is  not  obvious  from  the  very  beginning.  However,  a  cer¬ 
tain  set  of  potential  coefficients,  called  Blanchet-Damour  moments,  defined  to 
first  post-Newtonian  order  has  especially  attractive  features.  Moreover,  by  using 
such  Blanchet-Damour  potential  coefficients  we  get  the  simplest  possible  form 
of  the  multipole  expansion  of  the  post-Newtonian  potentials  (these  expansions 
take  an  almost  Newtonian  form).  Basically  two  sets  of  BD-moments  occur  in 
the  formalism:  mass-multipole  moments  and  (mass)  current  multipole  moments. 
Expressed  in  terms  of  (symmetric  and  trace-free)  Cartesian  tensors  they  are  de¬ 
noted  by  Ml  and  Sl-  Here  L  stands  for  a  Cartesian  multi-index,  L  =  i\  . .  ,i{ 
and  each  index  i  runs  over  the  three  spatial  indices.  The  set  Ml  is  equivalent  to 
a  set  of  potential  coefficients  C\m  and  S\m  that  appear  in  a  spherical  harmonic 
expansion  of  the  potentials.  The  first  spin-moment  of  a  body  agrees  with  its  spin 
or  total  angular  momentum.  The  multipole  expansion  of  We  and  WE  reads  (a 
dot  indicates  the  time  derivative): 
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where 


a  =  r\'-t±L*±±v  a  J_ 

^  (/  +  !)!  2/  +  3Pi  L|X|’ 


PL  =  /  S“  Xai  d3X, 

Jv 

and 

s(r,x)  =  r00  +  -4rss,  s“(r,x)=  -r0a, 

c2  c 

=  7*^(1,  ad)  are  the  components  of  the  energy-momentum  tensor  in  the 

GRS. 

The  function  A  does  not  enter  the  post-Newtonian  equations  of  motion.  The 
latter  contains  only  the  BD  multipole  moments  Ml  and  Sl-  The  only  place 
where  the  function  A  should  be  accounted  for  is  the  transformation  between  the 
proper  time  of  an  observer  and  the  coordinate  time  of  the  GRS.  First  estimates 
indicate  that  the  A-terms  will  be  negligible  in  the  foreseeable  future.  For  these 
reasons  the  gauge  function  A  will  not  be  mentioned  in  the  recommendations. 

A  spherical  harmonic  expansion  of  We  reads  ( R  =  |X|) 
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Here  the  relativistic  time  derivative  terms  are  not  expected  to  play  some  inter¬ 
esting  role  in  the  near  future;  they  will  not  be  mentioned  in  the  corresponding 
recommendation.  The  gravitomagnetic  vector  potential  of  the  Earth,  WB,  is 
dominated  by  the  Earth’s  spin- vector  Se  (total  angular  momentum),  i.e.,  to  a 
good  approximation 


Wf(T,  X)  ~ 


G(Xx  SE)“ 

~2 


-A, 


Note  that  this  representation  of  the  geocentric  metric  correctly  yields  the  Schwarz- 
schild  and  Lense-Thirring  accelerations  in  satellite  motion  to  first  PN-order  as 
recommended  e.g.,  by  the  IERS  Conventions  2000. 


3.5.  The  barycentric  metric  in  the  mass-monopole  approximation 

For  many  applications  it  is  sufficient  to  keep  the  mass-monopoles  of  the  various 
bodies  only,  i.e.  to  put 


Ml  =  0  for  /  >  1,  Sl  =  0  for  /  >  1 

for  all  bodies  and  to  keep  the  masses  Ma  =  Ma  only.  Furthermore,  we  will 
assume  all  moments  Vl  to  vanish.  From  the  transformation  rules  for  the  metric 
potentials  one  derives  the  metric  in  the  barycentric  coordinate  system  in  the 
form 


with 


9  oo  — 

9oi  =  — 


1  +  4W0 (t,  x)  -  4(wo(t,  x)  +  A (/,  x)), 


(19) 


where 
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A(t,x)  =  ^  ab(*,x), 
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Futhermore,  in  our  approximation 

wt(t,x)  =  J2^—^vB-  (23) 

B  TB 

Note,  that  the  post-Newtonian  Einstein-Infeld-Hoffmann  equations  of  motion 
for  a  system  of  mass-monopoles  that  form  the  basis  of  modern  solar  system 
ephemerides  can  be  derived  from  that  form  of  the  barycentric  metric.  Thus, 
the  barycentric  mass-monopole  metric  given  above  is  already  in  use  for  the 
description  of  solar  system  dynamics. 
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